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deterministic modeling of terrorism 
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Abstract 


A deterministic mathematical model of terrorism with government inter- 
vention was constructed from five compartments and subdivided into two 
core and non-core groups. A non-core group is a general group G(t), while 
the core group is susceptible S(t), moderate I(t), terrorism T(t), and re- 
covered R(t). The Elzaki transform method with differential transform 
to handle the nonlinear terms is employed to solve the model. The results 
show that government intervention on susceptible groups proved to be 90% 
effective in reducing terrorist threats since the group of susceptible in the 
population appears to be at risk of adopting the ideology through different 
means of contact. Also, due to the government intervention, the moderate 
group reduces gradually in time. 


AMS subject classifications (2020): 93A30, 65C20, 39A39. 
Keywords: Modeling; Terrorism; Elzaki transform; Differential transform 


and Intervention. 


1 Introduction 


Terrorism can be seen as extra-normal violence or a threat to intimidate 
innocent civilians through imposing their ideology, religiosity, and political 
objectives in a population. Global terrorism index 2020 defined terrorism 
as a threatened, or actual use of illegal force and violence by a non-state 
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actor to attain political, economic, religious, or social goals through fear, op- 
pression or intimidation. This definition deduces that terrorism is not only 
the physical act of an attack but also the mental blow it had on human- 
ity for many years; see [18]. Government intervention strategies consist of 
law enforcement. When violent extremists are investigated, prosecuted, and 
imprisoned, rehabilitation involves psychological counseling and correcting 
of their extremist views through religious education followed by vocational 
training. The military strategy consists of when violent extremists are killed 
or captured. Counter-terrorism, however, agrees that these strategies cannot 
reduce the threat of terrorism solitarily in [26, 8]. 


Similarly, the authors of [24, 25] stated that governments intervention 
could use an additional set of strategies collectively known as countering vio- 
lent extremism and classified it into three programs: 1) Prevention program 
which seeks to prevent the terrorist process from occurring, 2) disengagement 
program, which attempts to stop or control terrorist from occurring, and 3) 
terrorist recovery program, which attempts to change an individual extremist 
belief to re-integrate them into the society. These programs frequently aim 
to target convicted terrorist foot soldiers. Martcheva [22] discussed a mathe- 
matical model as a system using mathematical tools and language. Castillo- 
Chavez and Song [10] first developed transmission dynamics and spread of 
terrorism via mathematical modeling, in which they incorporated the spread 
of radical ideologies, fanaticism, recruitment, and terrorist activities. The 
models of terrorism were first developed in [10]. A mathematical model of 
the dynamic behavior of terrorism was developed in [1], and the authors con- 
trolled to check the recruitment pool. The model is developed to control the 
spread of terrorist ideologies in society. Furthermore, Santoprete [26] built 
a model of countering violent extremism (CVE) for prevention programs. It 
seeks to stop the radicalization process from occurring. According to [13], 
it was stated that Elzaki strictly connected with the Laplace transform has 
some advantages such as Elzaki over Laplace converges faster to the exact 
solution. 


Moreover, Akinola, Akinpelu, and Oladejo [6] applied the Elzaki decom- 
position method for solving the epidemic model. The obtained result agreed 
with the Adomian decomposition method, homotopy perturbation method, 
variational iteration method, and differential transform method (DTM). Sule- 
man et al. [31] used the Elzaki projected differential transform method to 
investigate the numerical solution of the fractional-order of the HIV model. 
According to [14], the Elzaki transform method (ETM) is a new integral 
transform, which is based on the Fourier series to handle only linear dif- 
ferential equations. Also, in [15, 23, 17, 16, 12, 9, 33], the Elzaki method 
transforms the fundamental properties and their applications to differential 
equations for the exact solution and analyzes boundary value problems via 
the Elzaki transform. The ETM was introduced to facilitate the process of 
solving ordinary and partial differential equations in the time domain. This 
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method is incapable of handling nonlinear equations, because of the difficul- 
ties that are caused by the nonlinear terms in [15, 23]. 

Furthermore, Kalavathi, Kohila, and Upadhyaya [21] presented the de- 
generate Elzaki transform to examine some of its properties and relations. 
Also, they investigated a scale preserving theorem for the degenerate Elzaki 
transform and theoretical dual transform to the degenerate Laplace trans- 
form. Aggarwal, Mishra, and Kumar [3] showed a comparative study of two 
new integral transforms Mohand and Elzaki transforms, respectively. Also, 
Anjum et al. [7] applied the ETM on the numerical solution for nonlinear 
oscillators. Bhadane and Ghadle [9] solved a system of linear differential 
equations using the ETM. Although, in [27], it was proposed a numerical 
solution for a third-order ordinary differential equation, using the DTM and 
ETM, then compared the results to see which method converges quickly to 
the exact solution. Some new applications of the ETM were defined in [29] 
for the solution to linear Volterra type integral equations. The results ob- 
tained were compared with solutions of the Adomian decomposition method, 
variational iteration method, method of successive approximation, Galerkin 
method, and Laplace transform method. The authors developed dualities 
in [2, 4] between the Elzaki transform and some useful integral transforms 
and then developed equations with decay problems. Their results showed 
that the Elzaki transform is in complete agreement with some useful meth- 
ods. Therefore, in [31, 20, 35, 28, 32], it was worked a new technique of the 
ETM for solving differential equations to an exact solution. Their obtained 
results showed that Elzaki converges to the exact solution faster than any 
other method. The Ezaki transform is a more powerful tool than the Laplace 
transform, which can still serve as an auxiliary method to the Laplace trans- 
form. 

This study aims to modify the existing models in [1, 26, 24] by incorporat- 
ing terrorist recovered groups and government intervention. Then we apply 
the ETM to solve the model equations together with DTM, to linearize the 
nonlinear term in the equations for numerical/exact solutions and compare 
the result of the methods on a graph. 


2 Materials and method 


Mathematical modeling of terrorism with government intervention is built, 
and then we apply numerical methods called ETM and DTM, respectively, 
to solve the model equations. 


2.1 Basic functions of ETM 


A set of function f(t) of exponential order is defined by A as follows: 
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A= {4 :4M,ky, ke > 0,|f(t)| < Mexp( a ifte (-1)10,00)} , (1) 
where M is a real number and k, and ko are integer. Then the Elzaki 


transform denoted by the operator EF as the new integral equation is defined 
by 


E|f(t)] w=u f re Jexp (= *) at, t>0,ky <u<ke,0O<t<o. 
0 
(2) 
Theorem 1. Let T(u) be the Elzaki transform. Then f(t) is a function 
taking the Elzaki f(t) as 
B[f'(®] = 7 -us. 
B[s"()] = 2 - s0 )— ¥F'0) 


E [f(a] = Pu) -¥ u2- aie 


Proof. Suppose that E [f =u i f'(t) exp (—+) dt. Using integration by 


part gives F [f’(t)] = Fu) ) — uf (0) mee Elf" (t)] = a — f(0)—uf’(0). The 
last expression of the theorem can be proved by mathematical induction. 


Therefore, Elzaki, Elzaki, and Elnour [17] stated the properties of Elzaki 
transform of a function below: 

Let T,(u) and T2(u) be the Elzaki transform of F(t) and F9(t), re- 
spectively. Then the Elzaki transform of [aF\(t) + bF2(t)] is given by 
[aT;(u) + bT2(u)], where a and 6 are arbitrary constants, and the inverse 
of Elzaki transform can be defined as f(t) = E~'(T(u)). 


The Elzaki transform of some functions is useful when applying on differ- 
ential equations as follows: 
(1) If f(t) =1, then E(1) =u f exp (—4) dt = u?. 
0 
(2) If f(t) =t, then E(t =u J texp (—4) dt =n’. 
(3) If f(t) =¢?, then E(t?) =u f @? exp (—£) dt = lu’. 
0 
(4) If f(t) =t",n EN, then E(t”) =u f t” exp (—4) dt = nlu”*?. 
0 
(5) If f(t) = exp(at), then E(exp(at)) = u f exp(at) exp (—+) dt = a 
0 
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2 


(6) If f(t) = cos(at), then E(cosat) = u f cosat exp (—+) dt = ——* 
0 


1-a2u2° 


1ta2u2* 


(7) If f(t) =sin(at), then E(sinat) = u f sinat exp (—+) dt = —2au" 
0 


(8) If f(t) = tcos(at), then E(t cos at) = u f tcosat exp (—+) dt = oe. 
0 


2.2 Advantages of Elzaki transform in comparison with 
Laplace transform 


In this subsection, the advantages of ETM, in comparison with the Laplace 
transform method are as follows: 


(a) The ETM has the capability of combining two powerful methods at the 
same time, compared to the Laplace transform method. 


(b) The Ezaki transform method is a powerful tool compared to the Laplace 
transform, which can still serve as an auxiliary method to the Laplace 
transform method. 


(c) The ETM converges to exact solutions faster than the Laplace trans- 
form method. 


(d) The ETM is a new integral transform that does not require perturba- 
tion. 


2.3 Differential transform method (DTM) 


There are many methods for solving differential equations. One of them is the 
Taylor series. The new form of the Taylor series called DTM was proposed 
by [34], and they applied it to solve mathematical problems in an electrical 
circuit. The idea of DTM is to determine the root of the Taylor series of a 
function by solving the induced recursive equation. In [19], it was developed 
a mathematical model for the transmission dynamics of the Syphilis disease 
and employed the DTM to handle the nonlinear systems. The DTM is a semi- 
analytic technique that depends on the Taylor series. The Taylor polynomial 
of degree n is defined by 


Pr(z) = > ALF (Ole — 0). (3) 
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Suppose that the function f has (n +1) continuous derivative on the interval 
(c—r,c+r),for all r > 0, and that lim, ... Rn(x) = 0, where R,,(x) is the 
error between P,,(x) and the function f(a). Then the Taylor series expanded 
about « = c converges to f(x); that is, 


fa) = ¥ AUOOe- oN, forall 2 € (c—r,c+r). (4) 


Differential transform of the function f(a) for the kth derivative is given in 
[5, 19] as 


F(k) = At {2 Jexa0; (5) 


dak 
where f(z) is the original function and F(k) is the transformed function. 
Using the ideas of [5, 19], the inverse differential transform of F'(k) is defined 
as 


Fla) = ¥. (a ~ 20)* Fb) (6) 
Substituting (7) into (8) yields 
f@)= S(e- m4 [2]. (7) 
k=0 xL=Xo 


The basic idea of DTM is considered in Table 1. 


Table 1: The fundamental operations of DTM 


Original Function Transform Function 
y(x) = u(x) + m(a) Y(k) = U(k) + M(k) 
y(a) = am(z) Y(k) = aM(k) 
y(«) = ae Y(k) = (k+1)U(k +1) 
y(z) = SE) Y (k) = (k+ 1)(k + 2)U(k +2) 
ye) = EVR) = (H+ DHF 2K + MUCK +n) 
y@)= 1 ¥ (R= oth) 
y(a) = 2 Y(k) = d(k — I) 
y(“z) =a™ at) = 0(k —m) 
y(x) = g(x)h(x) Y(k) = IGE —m) 
y(at) = exp (Ae) ¥(k) = i 
y(e)=(+2)™ Y(k) = eee 
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3 Formulation of the model 


The flow diagram of terrorist network within a host community can be de- 
duced in Figure 1. 


UR (4+ 66)T 


Figure 1: Flow diagram of terrorist with government intervention 


The model of terrorism with government intervention corresponding to 
the flow diagram in Figure 1 is given by the deterministic system of nonlinear 
differential equations as 


( 
tt Bs ( ay on S(E)I(t) — (w+ w+ 8)I(6) — Ba(ANT(E) (8) 
© = BAI) — (n+ OT) — ATU 
a = \T(t) + wl (t) — wR(t) — oR(t) 


The model focuses on a population of interest T(t), which is divided into 
noncore group G(t) and core group C(t) in [11]. The noncore group Git) 
is the general population, which the individuals in the group are at risk 
of adopting the ideology. The noncore group is usually the source of the 
recruitment pool, where the entire core group builds their members from. 
The core group C(t) consists of susceptible group S(t), the moderate group 
I(t), Terrorist group T(t), and recovered group R(t). The susceptible group 
S(t) is individuals who have not yet been converted into the ideology but 
have begun to walk around and fall into terrorist beliefs. Moderate group I(t) 
includes individuals who have been transformed with beliefs to reinforce their 
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terrorist views (Training stage). The terrorist groups T(t) are individuals 
who have completely converted and have received training to become either 
leaders or foot soldiers of the terrorist respectively. The recovered group R(t) 
is individuals who are neutralized due to counter-terrorist /NGOs activities 
applying to each group. 


Table 2: Description of parameters of terrorism model equation (10) 


Parameters Description 
A Recruitment per-capital (birth/migration rate) 
0) Death due to suicide attack or arrest by counter-terror activities 
Bi Rate at which terrorists recruit their members 


from general to moderate group 

Bo Contact rate between susceptible and moderate groups 
Lb Natural death rates 

Bs Contact rate between moderate and terrorist groups 
A 

x 

N 

w 


Probability the government measure intervention on each groups 
Immigration rate 
Rate at which individuals recover from the ideology 
Probability that the moderate terrorists move to R(t) 
due to government intervention 
Rate at which terrorist individuals progress 
from R(t) to G(t) after being certified by the counter-terrorist /NGO’s 


Table 3: Description of variables of terrorism model equation (10) 


Variables Description 
G(t) Number of general group at time t 
S(t) Number of partial-terrorist at time t 
I(t) Number of moderate terrorist (indoctrinization) at time t 
T(t) Number of terrorist at time t 
R(t) Number of terrorist recovered at time t 
We have 


C(t) = S(t) + 1(t) + T(t) + Rit). (9) 

The total population is 
N(t) = G(t) + C(t). (10) 
Then (1, 2,83 are parameters, which define the recruitment pool of each 
subgroup with each other. Similarly, 6) (wis) G(t)C(t) is the transition 
rate from G(t)-S(t). The moderate group has more chances of recruiting, 


since the group builds its members from S(t) and transfers them into the 
terrorist group individuals. There will accept their duties and execute the 
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mission given to them (to become complete terrorist at T(t)), in [11]. The 
dynamics between S(t)-I(t) and the transfer rate from I(t)-T'(t) within the 
core are [5 (cts) S(t)I(t) and 63AI(t), respectively. 

The intervention rate A is bounded in an interval 0 < A < 1, which is 
the measure of intervention on each compartment. The intervention rate A 
is the military/NGO’s activities, which aimed at reducing the recruitment 
pool/threats, respectively. Due to the intervention, some individuals will 
move to a recovered group. After some time, individuals who received either 
jail time or certified repentance and others will fall back to G(t). It is assumed 
that recovery from terrorism is not permanent, since it deals with ideology, 
marginalization and poverty. Moreover, A is the birth/migration rate, y is 
the immigration rate, is the natural death from each compartment, 6 is 
terrorist induced death rate due to military intervention or by suicide, while 
A is the rate at which individuals recover from the ideology due to intervention 
and ¢ is the rate at which terrorist individuals progress from R(t)-G(t) after 
being certified by the counter-terrorist /NGO’s. 


3.1 Assumptions of the model 


From the model formulated, it is assumed that A is the intervention rate, 
while C(t) is equal to 1, (C(t) = 1). Similarly, it is assumed that N(t) = 1. 
The recovered rate of individual terrorists in this model is not permanent 
since their behavior is dynamics. This means individual terrorists can change 
after being integrated into society can be re-engaged. Also, A is assumed 
to be a probability that not all individuals terrorists the intervention will 
capture (some will hide). The rate at which an individual leaves the terrorist 
group to enter the recovered group is A, where @ is the fraction at which 
terrorists successfully move to the general group due to the intervention by 
the government. Susceptible, moderates and terrorists are assumed to have 
frequent contacts with the general group. 

Finally, the age limit of an individual ranges from five years and above 
while everybody in the population has the same average natural death rate 


[Lb 


3.2 Solution of the model using ETM and DTM 


Consider the nonlinear differential equation corresponding to the model dia- 


gram as 
dG 
Hp TAT Yt+ OR (Ht BANG, (11) 
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dS 

ae Bi AG — pS — (A) ST, 
dl 
tea B2(A) SI — (w+ ut+d+ B3A)l, 


aT 
p= BAL (Ut 8+ NT, 
dR 
dt 


So, using (14)-(15) gives the following equations: 


= AT +wl— (w+ @)R. 


dG 
dt 
dS 


ae = [3G — pS — LIS, 


=1,+¢R— L2G, 


dl 
—=14IS—LsI 
dt 4 545 


dT 
dt 
dR 


Gy TAT + wl — LR, 


ae FS enw Fret 


where 


[,=A-y, 
Lg =p+ PiA, 
L3 = Bi A, 
[4 = BoA, 

Ls =wt+pt+ot 3A, 
Lg = B34, 
L7=pt+o+ A, 
Lg =p+¢. 


Applying the ETM on both sides of equations (16)—(20) gives 


E (=) = E{L,+¢R- L2G}, 
E (F) = E{L3G — uS — Lal S}, 


E (=) = E{L4IS — LsI}, 


dT 


E (=) = Pi Lel iat} 
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(12) 
(13) 
(14) 


(15) 


(16) 
(17) 
(18) 
(19) 


(20) 


(21) 


(22) 


(23) 


(24) 
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E (F) = Hie furs teks (25) 


Suppose that E [G(t)] = T\(u), that E [S(t)] = To(u), that FE [I(t)] = T3(w), 
that F [T(t)] = Ti(u), that E[R(t)] = T5(u) and that Ly in (21) is a con- 
stant. Then using the ETM properties of constant on L, yields E(L,) = 
u fy La exp(=)dt = w{ Jim (ZL; exp(=*) +ul1)} = u*L; as E(Li) = 
uy. 


Applying the DTM on the nonlinear part in (22) and (23) implies 


£(B) = Bla — pS ~ Lin} (26) 
dl 
E (=) ny) 5. eee Fe aS (27) 
where A, = > S(m)I(k — m), which is defined by the DTM and the poly- 
nomial as 
Ao = Solo, 
Ay = Soli + SiIo, 
Az = Sola + Sif, + Solo, (28) 
As = Sol3 + SyIg + Solo, 
Subsequently 
[A — co] = bale] + ofa (u)]  LalTiCu)) (29) 
J) — uso] = Zaft] —nlTa(w)] BULaAnl (80) 
[EOD — ra(0)] = Leltalw)] — Lala) (32) 
[A — uro(0)] = ATa(w)] + tTa(u)]~ Eafe), 83) 


with the initial conditions 


G(0) = Go(t) = Go = 100, 
S(0) = Sy(t) = So = 50, 
T(0) = To(t) = To = 25, 
R(0) = Ro(t) = Ro = 10, 
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By simplifying and collecting the same terms of (29)—(33), we have 


F 4: 1s] Pili) OTR (G) =A by eas, (35) 
[5 + a] To(u) ~ LoTi(w) = Blan] + us, (36) 
F 1s] Pi) SBA Palys (37) 
E 4 1] Ti(u) — Le T3(u) = uTo, (38) 
E + ts| T5(u) — wT3(u) — AT4(u) = wRo. (39) 


Solve equations (35)—(39), for Ti(u), To(u), T3(u), Ta(u), and Ts(u). Then 
truncate A, at So, Io, using (28) and applying the ETM. Subsequently, the 
DTM is employed to handle the nonlinear terms in (22) and (23), and we 
solved the system with the help of Maple 17 software (Intel (R) Celeron (R) 
CPU@ 1.10GHz, (RAM) 4.00GB), which gives the following results: 

1 
(14+ ubg)(uL2 +1)(1 + ub7)(uls +1) 
+? LgliLsL7 + (((Li + Rod + GoLs)Ls + LiLg + dwlo)L7 
+($XTb + LiLg)Ls + ¢AL61o)u? 
((Rod + L5sGo + Gols + Li) L7 + (L1 + Rod + GoLs)Ls 
+ (ATo + wlo)o + LiLs)u? 
(Li + GoL7 + Rod + LsGo + GoLs)u + Go)u?) (40) 


((uALg + woL7) + woBuLa(An) 


Ty (u) => 


TG) CO ihs E al ELS DOT AC a aie Pag be 
+D3¢w + LoLgLs)u? + ((Lg + Ls + Lo)L7 + (Lg + La) Ls + L2Lg)u? 
+(Ds + Do + L7 + Lg)u)EL4(Ay) 
u(w? L3LgLiLsL7 + (((Li + Rod + Gols) L3 + LgL2So)Ls 
+ [3(L1Lg + dwly))L7 
+L3((ATp + LiLs)Ls + oAL6Io))u* + (((L3Go + So(Ls + L2))Ls 
+(Li + Rod + Gols) L3 + LgL2So)L7 
+((Li + Rog + GoLg)L3 + LgL2So)Ls + L3(LyLg + (AT + wlo)¢))u° 
+((SoL5 + L3Go + So(Lg + Le))L7 
(L3Go + So(Lg + L2))L5 + (Li + Rod + GoLg)L3 + LgL2So)u? 
+(SoL7 + Sols + L3Go + So(Lg + L2))u + So))u)/ 
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(1+ uLg)(uL, + 1)(1+ uL7)(uLs + 1)) 


(41) 
_ (EL4(An) + ulo)u 
LE) meer (42) 
_ u?(EL4(An)Le + (L5To + Lelo)u + To) 
TAOS (1 + uL7)(uLs +1) ~ 
Ts(u) = (i+ yi 11+ ule) (u?(ALe + wL7)u + w)EL4(An) 

+ ((L5Ro + wl) L7 + ML5To + Lelo))u? 

i (RoL7 + XTo + Is Ro + wIo)u + Ro)). (44) 


Thus the initial condition and parameters are 

G(0) = 100, S(0) = 50, 7(0) = 15, T(0) = 25, R(O) = 10, A = 600, d = 
0,0022, w = 0.015, A = 0.008, A = 0.7, 61 = 0.0035, G2 = 0.000005, 63 = 
0.00045, 6 = 0.0036, 4 = 0.000034247 and y = 0.7. 


Applying the conditions above on (40)—(44) gives 


T,(u) =100u? + 599.0035754u3 — 1.90652029u* + 0.00605328738u° 
— 0.0000189756893u° +--+, 
T2(u) =50u? + 0.4708656101u> + 1.886410919u* — 0.00606738403u° 
+ 0.0000192353168u° + --- , 
T3(u) =14u? — 0.2800526593u? + 0.005228652322u* 
— 0.00009762023033u° + 0.000001822593812u° + --- , 
Ti(u) =25u? — 0.2847812124u3 + 0.003197551982u4 
— 0.00003499897343u° + 3.652596885107- “ue +--- , 
T5(u) =10u? + 0.3886565195u? — 0.007974537753u" 

+ 0.0001361224536u° — 0.000002327546774u® +--- . 


Take the inverse of the Elzaki transform to solve the induced recursive equa- 
tion of G(t) = E-1(Ti(u)), $(t) = E-*(Ta(u)), I(t) = E-*(Ta(u)), T() = 
E~1(Ty(u)), and R(t) = E~+(T5(u)). 

Then using Maple 17 software to solve the above expression, the closed form 
of the solution when k = 4 can be written as 


G(t) =100 + 599.0035t — 0.95326t? + 0.001008¢° + 7.9065 x 10-7#4+---, 
+ 0.0839¢7 — 0.000009706t? + 5.667 x 107 #4 + --- 
=15 — 0.28005t + 0.002614? — 0.00001626¢° + 7.594 x 10-S¢4+.--- , 


’ 


~—~"_—_—_~" ~—~—"_— ~~" 
II 
oO 
oO 
+ 
S 
j=) 
I) 
aD 
iw) 
(8) 
is 
| 
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R(t) =10 + 0.3886¢ — 0.003982? + 0.0000226¢? — 9.698 x 10-8#4 +--- , 


Furthermore, the DTM is also employed to solve (16)—(20) numerically 
and compare the result on the graph with that of Elzaki transform. The 
initial condition and the parameter values are 
G(0) = 100, S(0) = 50, 7(0) = 15, T(0) = 25, R(O) = 10, A = 600, ¢ 
0,0022, w = 0.015, A = 0.008, A = 0.7, 61 = 0.0035, G2 = 0.000005, 63 = 
0.00045, 6 = 0.0036, 4 = 0.000034247 and y = 0.7. 

Also, using the transformation in Table 1, for (16)—(20), gives 


I 


_ Li(k) + 6R (k) = L2G (h) 


G(k +1) arm] ; (45) 
I4 > S'(m) I (k — m) + L3G (k) — pS (k) 
Sh). = — et (46) 
ba > S(m)I(k—m) — Ls (k) 
I(k+1) =—™= a (47) 
r(k-+1) =O (48) 
Hove A) st nS (49) 


Applying the initial condition in (45) to (49) and using Maple 17 software, 
we have 


4 
G(t)= 5° G(e)t* 


m=0 


= 100 + 599.2905t — 0.09373t? + 0.0000069105t? + 1.1743 x 1078t4 +--- , 


4 
S(t) = >> Skye 


m=0 
= 50 + 0.02658t + 0.08389? — 0.000009522¢ — 1.8694 x 1079#4 +--+. , 
4 


rO= yo Tee 


m=0 


= 15 — 0.27975t + 0.002608t? — 0.00001605¢° + 7.24869 x 107-844 +--- , 


4 
Posy Pee 


m=0 


= 25 — 0.29031t + 0.0016837t7 — 0.000006498¢° + 1.87567 x 107844 +... r 
4 
R(t)= > R(k)t* 


m=0 


= 10 + 0.3886¢ — 0.003965t? + 0.000022338¢° — 9.3485 x 1078444... . 
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4 Numerical simulation 


The simulations were carried out using the following variables and parameters 
for initial conditions. The final time was t = 5 years. Computations were 
run in Maple 17 software (Intel (R) Celeron (R) CPU@ 1.10GHz, (RAM) 
4.00GB) for analyzing. Therefore, the variable and parameter values of source 
are 

G(t) = 100, S(t) = 50, I(t) = 15, T(t) = 25, R(t) = 10, A = 600 [26], 6 
= 0.0036 [24], 8: = 0.0035 [24], Bz = 0.000005 [24], 83 = 0.00045 [24], w = 
0.000034247 [24], A = 0.008 [26], A = 0.1, 0.7,0.9 varies, y = 0.7 [26], w 
0.015 and ¢ = 0.0022. 


General group 


Figure 2: Comparison graph between ETM and DTM on general group in time. 
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31 J —:-— 5M 


Figure 3: Comparison graph between ETM and DTM on susceptible in time. 
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Figure 4: Comparison graph between ETM and DTM on moderate group in time. 
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Figure 5: Comparison graph between ETM and DTM on terrorist group in time. 
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Figure 6: Comparison graph between ETM and DTM on recovered group in time. 
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Figure 7: Graph of susceptible group in time with different rate of government inter- 
vention A. 


5 Discussion of results 


The ETM was introduced to determine the approximate numerical solution 
of a system of ordinary differential equations. Maple 17 software was used 
for simulations. Government intervention on susceptible groups proved to 
be effective in reducing terrorist threats since the susceptible group in this 
research is seen to be at risk of adopting the ideology through contact. sim- 
ilarly, this can be deduced from Figure 3 with different rates of intervention 
A. 

The graph in Figure 2 for five years, showed that ETM and DTM agreed 
with each other without any significant difference in comparison. Similarly, 
this can be deduced from Figure 3 where the two methods also agreed on 
the comparison. Figure 4 shows the comparison graph between ETM and 
DTM on the moderate group against time, which is gradually reduced due to 
government intervention. From the graph, this can be seen when time dura- 
tion continues after 5 years, the methods for analytic solution may not agree 
with each other due to control on moderate terrorist since it is dynamics. 
Furthermore, a similar scenario can be deduced from Figure 5 on the com- 
parison between ETM and DTM. This can also be observed from the graph 
in Figure 6, where the effect of the comparison showed no significant differ- 
ence between ETM and DTM on the recovered group. Finally, Figures 2-3 
showed that both methods are more effective, more powerful, and accurate 
for solving a system of differential equations for an approximate solution, but 
Elzaki transform converges faster than DTM. 
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6 Conclusion 


A mathematical model of terrorism with government intervention was de- 
veloped to reduce the recruitment pools within the five compartments. The 
ETM and DTM were tested and solved the governing differential equations 
compared to the results on the graph. Government intervention on suscep- 
tible groups proved to be 90% effective as well as a moderate group. The 
susceptible and moderate group in this study appears to be exposed to adopt- 
ing the ideology through contact. The ETM was employed to solve the model 
with Maple 17 software and applied DTM to decompose the nonlinear term 
since some differential equations have difficulty finding analytic/exact solu- 
tions. Finally, a plan to incorporate other issues would be addressed in the 
future study. 
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